CASIMIR FORCE BETWEEN A SPHERE AND A PLANE: 
SPECTRAL REPRESENTATION FORMALISM 



CECILIA NOGUEZ AND C. E. ROMAN- VELAZQUEZ 
Instituto de Fisica, Universidad Nacional Autonoma de Mexico, 
do Postal 20-364, Mexico D.F. 01000, Mexico 



We develop a spectral representation formalism to calculate the Casimir force in 
the non-retarded limit, between a spherical particle and a substrate, both with 
arbitrary local dielectric properties. This spectral formalism allows one to do a 
systematic study of the force as a function of the geometrical variables separately 
from the dielectric properties. We found that the force does not follow a simple 
power-law as a function of the separation between the sphere and substrate. As 
a consequence, the non-retarded Casimir force is enhanced by several orders of 
magnitude as the sphere approaches the substrate, while at large separations the 
dipolar term dominates the force. 



1 Introduction 

Recent experiments 1-4 to measure the Casimir force have been done choos- 
ing a spherical surface and a plate, instead of the plate-plate configuration 
originally proposed by Casimir™. The interpretation of the Casimir force in 
these precision experiments has relied in the Proximity Theorem, developed by 
Derjaguin and collaborator iM The Proximity Theorem estimates the Casimir 
force per unit area between two curved surfaces of radii R\ and R2, in terms 
of the Casimir energy per unit area between parallel planes^, V(z), where z 
is the distance between plates. Assuming that the Casimir force on a small 
area of one curved surface is due to locally "flat" portions on the other curved 
surface, Derjaguin found that, 

= 2tt (-fi^A V(z). 



R± + i?2 , 

In the limit, when Ri = R, and R2 — > 00, the problem reduces to the case of 
a sphere of radius R and a flat plane, that yields to Fpt(z) — 2ttRV(z). The 
force obtained is a power-law function of the separation, Fpt(z) cx , where 
at "large" distances (3 = 3 while at short distances (3 = 2. This theorem is 
assumed to hold only when z <§; R\,Ri, however, it is not clear up to what 
limit this approach is valid. 

The origin of dispersive forces between atoms and macroscopic bodies 
may be attributed to interactions between the induced charge distributions 
on them by quantum vacuum fluctuations^. The charge distribution can be 
represented, in the simplest approximation, by electric dipoles. This dipolc 
approximation was employed by London 121 to calculate the van der Waals 
energy V v w(z) between two identical polarizablc molecules by using pertur- 
bation theory in quantum mechanics. Considering that the molecules interact 
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through a quasi-static potential (or non-retarded limit), London found that 
Vy-w(z) oc — l/z^ 1 where (3 — 6 and z is the magnitude of the distance between 
molecules. In 1948, Casimir and Polder 121 studied the interaction of a neutral 
atom with a perfectly conducting plane, such that, the force is proportional 
to —l/z@ with p = 4 in the non-retarded limit El Furthermore, they found 
that the vW interaction could be attributed to the change of the electromag- 
netic zero-point energy of the classical proper electromagnetic modes of the 
system. The same behavior can be also obtained calculating the change of 
the zero-point energy between a induced dipole moment on a sphere of polar- 
izability a with its own image dipole in a plane. Contrary to what one may 
expect, the result for perfect conductors differs from the one obtained using 
the Proximity Theorem by a factor of 1/z 2 . Therefore, an exact calculation 
of the Casimir force between a sphere and a planar surface becomes essential. 

In this paper, we study the Casimir force in the non-retarded or quasi- 
static limit between a sphere and a plane, both with arbitrary dielectric prop- 
erties. We calculate the difference of the zero-point energy when the two 
bodies are at a distance z, and when they are at infinite, as 

£(*) = -£>*(*)-<,], (i) 

s 

where lu s (z) and w^, are the proper electromagnetic modes at z and z — ► oo, 
respectively. By including multipolar interactions to all orders, we calculate 
exactly the proper el ectr omagnetic modes of the system using a Spectral Rep- 
resentation formalisirffl "We show that this spectral representation formalism 
has the advantage that separates the contribution of the dielectric properties 
of the sphere from the contribution of its geometrical properties, allowing to 
perform a systematic study of the system as a function of the radius of the 
sphere and its dielectric properties. 

2 Dipole approximation within the non-retarded limit 

Let us consider a sphere of radius a and dielectric function e s (uj) which is at a 
minimum distance z from a semi- infinite plate of dielectric function e p (u>), as 
shown in Fig. 1. The sphere and plate are electrically neutral. The quantum 
vacuum fluctuations will induce a charge distribution on the sphere which also 
induces a charge distribution in the plate, such that, the Zm-th multipolar 
moment on the sphere is given by 

Qlm = Otlm [V/™ C + V^" b ] , (2) 

where VJ™ C is the field associated to the quantum vacuum fluctuations at the 
zero-point energy, VJ™ b is the induced field due to the presence of the plate, 
and aim is the Zm-th polarizability of the sphere. 

In the dipolar approximation, one can find using the image method, that 
the charge distribution on the plate is described by a dipole moment p p {oj) 
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Figure 1. Schematic model of the sphere-plane system. 



which is related with the induced dipole moment on the sphere by 
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PpH = i 



e p (u>) 



•Ps(w), 



(3) 



where M = (1,1,-1) is a diagonal matrix, and p s is the dipole moment on 
the sphere which, according with Eq. J2J , is given by 



ft(w)=a(w) E vac + T-p p {u) 



(4) 



Here, T = (3rr — r 2 I)/r 5 is a tensor that couples the interaction between 
the induced dipoles in the sphere and plate, 1 is the identity matrix, and 
r = (0, 0, 2(z + a)) is the vector from the center of the sphere to the center 



of the image charge on the plane, such that, 
Substituting Eq. J3J in Eq. (gj, we find 



T = (-l/r 3 ,-l/r 3 ,-2/r 3 ). 



E vac . (5) 



with / c = [1 — e p (oj)]/[l + e p (w)], and b = 1, 1 or 2. The eigenfrequencies 
of the sphere-substrate system are found from the above equation and must 
be independent of E vac . The normal frequencies can be obtained when the 
determinant of the matrix in the left-side of Eq. (JSJ) is equal to zero. Then, 
there are three modes (two are degenerated), with frequencies given by the 
solution of 

n 2 r 



1 



+ fc 



1 



[2(z + a) 



1 



= 0. 



(6) 



One should notice that the left-side of each term of the above equation is 
related only with the dielectric properties of the sphere, while the right-side 
is related with its geometrical properties only. 
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3 Spectral Representation of the sphere-plane configuration 

The polarizability of a sphere in the non-retarded limit is a /a 3 = a — [e s — 
l]/[e s + 2]. Now, let us define the variable u = [1 — e s ] _1 , then a(oj)~ l = 
—3u(uj) + 1, and substituting in Eq. ©, one obtains 

[- u Hl + H]- Ps H = V vac , (7) 

with V vac = ±a 3 E vac , and H = | [I — / c a 3 T • M]. The later is a dimensionless 
matrix that only depends on the geometry of the system through the ratio 
z/a. 

Let consider when f c is reaff^, then H is real and symmetric, in such 
case, we can always find a unitary transformation U that diagonalizes it, 
U -1 lHnU = n s , being n s the eigenvalues of H. Therefore, the proper modes 
depend on the separation z, and can be obtained from 

det [-u{u)I + H] = + n s (z)] = 0. (8) 

s 

Alternatively, the eigenfrequencies can be found from the poles of the Green's 
function of the system defined as <G(u) = [— u(w)l + H] _1 , whose ij-th element 
is given by 

G^u) =Y, Uis{UjS Tl ■ ( 9 ) 
^— ' u — n s (z) 

Once the proper modes are obtained, the interaction energy is calculated using 
Eq. (1), as the difference between the zero-point energy when the sphere is at 
a distance z from the substrate and when it is at z — > 00, then the Casimir 
force is calculated as F — —d£(z)/dz. 

One should notice that in the spectral representation formalism, the mate- 
rial properties of the sphere are contained in the spectral variable u, while the 
geometrical properties, like the radius of the sphere and the separation of the 
sphere to the substrate are in H. Furthermore, H is dimensionless and depends 
on the ratio z/R, its eigenvalues are independent of V vac . On the other hand, 
the dielectric properties of the substrate are in f c which can be a real function 
even for dispersive materials. Therefore, with the spectral representation one 
can study separately the contribution of the dielectric properties of the sphere 
from its geometrical properties, allowing for a systematic investigation of the 
system. 

4 Exact calculation within the non-retarded limit 

The ^m-th multipolar moment induced in the sphere on Eq. (2) can be written 
using a multipolar expansion-^, that yields to 

Qim = ~^ l) ai m [V^ + {-IT +l ' A*™' Q Vm ,], (iO) 

I' ,m' 
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where Q/'m' is the I'm'-ih induced multipolar moment in the substrate located 
at r = (2z + 2a, it, ip) from the center of the sphere, and A\™ is the matrix 
that couples the interaction between sphere and substrate^. The induced 
/'m'-th multipolar moment in the substrate is related with the multipolar 
charge distribution on the sphere like Qi'm' = {~ 1)' +m fcQi'm'- Substituting 
this in Eq. (|10fl . one finds 



E 

I'm' 



^SwSmm' , w'm' 



{21 + l)a v 



Qi'm' = VZ C - (11) 



If the sphe re is homogeneous its polarizabilities are independent of m, and 
are given bjU^ 

„,(, A - *[ £s M - 1] 21+1 _ n io n 2i+i n9 N 

<xi(v) — 77 , r , a - —a , (12) 

l[e s (u>) + lj + 1 n/ -w(w) 

where n; = Z/(2Z + 1). The poles of the right-side of Eq. (JT^J, u(cj; ) = 
n;o, yield the frequencies of the proper modes of the isolated sphere. Again, 
we have separated the material and geometrical properties of the sphere in 
Eq. C2J. 

Using Eq. (fT2"|) . we can rewrite Eq. © as 

^[-)i( u )V+^']v = ^ 

Q !m (Za 2Z+1 ) 1/2 
where /z = (Z,m), x M = ^ q2 ; + 7)1/2 ' *V = — V{^ c , with 



Hji' = nroV + / c (a '^ +1) (/n 1/2 <- (13) 

It was shown in Ref. 12 that, ^ is a symmetric matrix that depends on the 
distance between the center of the "image-charge" distribution and the center 
of the sphere as, l/[2(z + a)] l+l +1 . Then, H is dimensionless, symmetric, and 
depends only on the geometry of the system. As we showed, it is possible to 
find the Green's function, G, of the system, in terms of the unitary matrix that 
diagonalizes H. The poles of G, yield the proper frequencies of the system. 



5 Results and discussion 



To illustrate the procedure, we use the Drude model for the dielectric function 
of the sphere, e s (uj) = 1 — Wp/[ui(w + i/T)], with us p the plasma frequency and r 
the relaxation time. In this case, u(w)w 2 = uj(lu + i/r), such that, the proper 
modes are given by 

u a {z) = -i/2r + J (i/2r) 2 + u%n s {z) s» u p ^Jn s {z), (14) 
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with s = (/, m). In the left-side we have considered that 1 ^> (tuj p ) 1 . In this 
case, the fluctuations of the zero-point energy, according with Eq. (1), are 

£ ( z ) = -7T ^ W^i^ - V™io\ ■ ( 15 ) 

One should notice that it is possible to find the behavior of the energy as a 
function of z/a for any Drude's sphere independently of its plasma frequency, 
like £ = £/hujp. Furthermore, the Casimir force can be also studied indepen- 
dently of the plasma frequency and radius of the sphere, since £(z/ a) one can 
find the force like 

d£(z/a) d£(z/a) d(z/a) 
dz d(z/a) dz 

and then define the dimcnsionless variable Fa = Fa/htOp. In summary we 
have shown, for a Drude sphere with 1 3> (tcj p ) _1 , that given the dielec- 
tric properties of the substrate the behavior of the energy and force is quite 
general, showing the potentiality of the Spectral Representation formalism. 




Figure 2. (a) Energy E, and (b) force Fa as a function of z/a. 

In Fig. 2, we show the energy £ and force F a fo r a Drude sphere, of 
arbitrary plasma frequency and radius over sapphir^Sl as a function of z/a. 
We calculate the energy and force when multipolar interactions up to / and = 
2000 are taken into account, as well as up to dipolar, and up to quadrupolar 
interactions are considered. When we consider multipolar interactions greater 
than dipolar ones (1,1' > 1), it is observed that £ and Fa do not follow a 
simple power-law as a function of the separation. We find that £ and Fa 
are proportional to l/(z/a) /3 , where (3(z) is a positive integer which is also a 
function of the separation. As the separation decreases (z — * 0) the power 
increases (f3 —> oo). Then, the force can suddenly increases more than four 
orders of magnitude as compare with the dipolar approximation, when the 
sphere approaches the substrate. While at large distances (7a > z > 2a) the 
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force can be obtained exactly if up to quadrupolar interactions are considered. 
And for z > 7a, the interaction between the sphere and the substrate can be 
modeled using the dipolar approximation only, like in the Casimir and Polder 
model. Finally, the increment of the force at small separations could explain 
the physical origin of the large deviations observed in the deflection of atomic 
beams by metallic surfaces, as well as some instabilities detected in micro 
and nano devices. However, specific experiments and calculations have to be 
performed to prove the latter. 
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